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ABSTRACT 

Weak shear maps of the outer regions of clusters have been successfully used to map 
the distribution of mass at large radii from the cluster center. The typical smoothing 
lengths employed thus far preclude the systematic study of the effects of galactic-scale 
substructure on the measured weak lensing signal. In this paper, we present two methods 
to infer the possible existence and extent of dark halos around bright cluster galaxies 
by quantifying the 'local' weak lensing induced by them. The proposed methods are: 
direct radial averaging of the shear field in the vicinity of bright cluster members and a 
maximum-likelihood method to extract fiducial parameters characterizing galaxy halos. 
The correlations observed for early-type galaxies on the Fundamental Plane are used 
to derive the scaling laws with luminosity in the modelling of cluster galaxies. We 
demonstrate using simulations that these observed local weak-shear effects on galaxy 
scales within the cluster can be used to statistically constrain the mean mass-to-light 
ratio, and fiducial parameters like the halo size, velocity dispersion and hence mass of 
cluster galaxies. We compare the two methods and investigate their relative drawbacks 
and merits in the context of feasibility of application to HST cluster data, whereby 
we find that the prospects are promising for detection on stacking a minimum of 20 
WFPC2 deep cluster fields. 
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1 INTRODUCTION 

Clusters of galaxies are the most recently assembled struc- 
tures in the universe, and the degree of observed substruc- 
ture in a cluster is the result of the complex interplay 
between the underlying cosmological model (as has been 
demonstrate d by many groups including 



Bird (1993) 



Evrard 



et al. (1993| ) and [West fc Bothun (199C| )) and the physical 
processes by which clusters form and evolve. Many clus- 
ters have more than one dynamical component in th e ve- 
locity structure in addition to spatial subclustering (|CoL 
less & Dunn 1996 , [Kricssler, Beers, fc Odewahn 1995 , [Bird 



1993| , [West fc Bothun 199Cl| and [Fitchett 1988[ ). Substructure 
in the underlying cluster potential and specifically the sub- 
clumping of mass on smaller-scales (galactic scales) within 
the cluster can be directly mapped via lensing effects. 

The observed gravitational lensing of the faint back- 
ground population by clusters is increasingly becoming a 
promising probe of the detailed mass distribution within a 
cluster as well as on larger scales (super-cluster scales). We 
expect on theoretical grounds and do observe local weak 
shear effects around individual bright galaxies in clusters 
over and above the global shearing produced by the 'smooth' 
cluster potential. While there is ample evidence from lensing 
for the clumping of dark matter on different scales within the 



cluster, the spatial extent of dark halos of cluster galaxies 
are yet to be constrained. The issue is of crucial importance 
as it addresses the key question of whether the mass to light 
ratio of galaxies is a function of the environment, and if it 
is indeed significantly different in the high density regions 
like cluster cores as opposed to the field. Moreover, it is 
the physical processes that operate within clusters like ram- 
pressure stripping, merging and "harassment" that imply 
re-distribution of meiss on smaller scales and their efficiency 
can be directly probed using accurate lensing mass profiles. 

Constraining the fundamental parameters such as mass 
and halo size fro m lensing effects for field galaxies was at- 
tempted first by Tyson et al. (1984 ) using plate material, 
the quali ty of which precluded any signal dete ction. More 
recently, Brainerd, Blandford, fc Small (1996 ) used deep 
ground-based imaging and detected the galaxy-galaxy lens- 
ing signal and hence placed upper limits on the m ean mass 
of an average field galaxy. [Griffiths et al. (1996 ) used the 
Medium Deep Survey (MDS) and HST archival data in a 
similar manner to extract the polarization signal. Although 
the signal is unambiguously detected, it is weak, and no 
strong constraints can yet be put on the mean profile of 
field galaxies, but the prospects are promising for the near 
future. 

On the other hand no such analysis has been pursued 
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in dense regions like clusters, and very little is known about 
the lensing effect o f galaxy halos super posed on the lensing 



effect of a cluster. Kneib et al. (1996) have demonstrated 
the importance of galaxy-scale lenses in the mass modeling 
of the cluster A2218, where the effect of galaxy-scale com- 
ponents (with a mean mass to light ratio ~ 9 in the R-band) 
needs to be included in order to reproduce the observed mul- 
tiple images. Mass modeling of several other clusters has 
also required the input of smaller-scale mass components to 
consistently explain the multiple images as well as th e ge- 



ometry of the arcs, for instari ce , in the case of CL00 24 ( Kas- 
siola, I ^ovner, fc Fort (1993| ), ^mail et al. (1995a|) ), where 
the length of the three images of the cusp arc can only be 
explained if the two nearby bright galaxies contribute mass 
to the system. This strongly suggests that the dark mat- 
ter associated with individual galaxies is of consequence in 
accurately mapping the mass distribution, and needs to be 
understood better, particularly if clusters are to be used as 
gravitational telescopes to study background galaxies. 

The observed quantities in cluster lensing studies are 
the magnitudes and shapes of the background population 
in the field of the cluster. To reconstruct the cluster mass 
distribution there are many techniques currently available 
which allow the inversion of the distortion map into a rela- 
tive mass map or an absolute ma ss map if (i) multiple arcs 
are observed ( Kneib et al. 199(:) and or (ii) magnificatio n 
effects are included ( Broadhurst, Taylor, fc Peacock 1995 ) 



Recent theoreti cal work (iKaiser fc Squires 1993 , Kaiser 1995 



Schneider 1995 , Schneider & Seitz 1995 and Squires & Kaiser 
1995|) H as focused on developing various algorithms to re- 
cover the mass distribution on scales larger than 20-30 arc- 
sec, which is roughly the smoothing scale employed (cor- 
responding to ~ 100 kpc at a redshift of 2; ~ 0.2). These 
methods involve locally averaging the shear field produced 
by the lensing mass, and cannot be used to probe galaxy- 
scale perturbations to the shear field. 

Our aim in this paper is to understand and determine 
the parameters that characterize galaxy-scale perturbations 
within a cluster. In order to do so, we delineate 2 regimes: 
(i) the 'strong' regime where the local surface density is close 
to critical (k ~ 1, where k is the ratio of the local surface 
density to the critical surface density) and (ii) the 'weak' 
regime where the local surface density is small (k < 1). The 
'strong' regime corresponds to the cores of clusters, and in 
general involves only a small fraction (typically 5-20) of the 
cluster galaxies whereas the 'weak' regime encompasses a 
larger fraction (~ 50-200). We are restricting our treatment 
to early-type (E & SO's) bright cluster galaxies throughout. 

We compare in this analysis the relative merits of our 2 
proposed methods; a direct method to extract the strength 
of the averaged local shear field in the vicinity of bright clus- 
ter galaxies by subtracting the mean large-scale shear field, 
and a statistical maximum likelihood method. The former 
method affords us a physical understanding, helps estab- 
lish the importance and the role of the various relevant pa- 
rameters and yields a mean mass-to-light ratio; the latter 
permits taking the strong lensing regime and the elliptic- 
ity of the mass of galaxy halos into account correctly. Both 
approaches are investigated in detail in this paper using nu- 
merical simulations. 

The outline of the paper is as follows. In Section 2, we 
present the formalism that takes into account the effect of 



individual galaxy-scale perturbations to the global cluster 
potential. In Section 3, the direct method to recover these 
small-scale distortions is outlined and in Section 4 we present 
the results of the application of these techniques to a simu- 
lated cluster with substructure. In Section 5, we examine the 
constraints that can be obtained on the parameter space of 
models via the proposed maximum-likelihood method. We 
also explore the feasibility criteria for application to cluster 
data given the typical uncertainties. The conclusions of this 
study and the prospects for application to real data and fu- 
ture work are discussed in Section 6. Throughout this paper, 
we have assumed Ho = 50kms~^Mpc~^, Q = 1 and A = 0. 



2 GALAXY-SCALE LENSING DISTORTIONS 
IN CLUSTERS 

2.1 Analysis of the local distortions 

The mass distribution in a cluster of galaxies can be mod- 
eled as the linear sum of a global smooth potential (on scales 
larger than 20 arcsec) and perturbing mass distributions 
which can then be associated with individual galaxies (with 
a scale length less than 20 arcsec). Formally we write the 
global potential as: 



(1) 



where <j!>c is the smooth potential of the cluster and (pp^ are 
the potentials of the perturbers (galaxy halos). Henceforth, 
the use of the subscripts c and p refer to quantities com- 
puted for the cluster scale component and the perturbers 
respectively. The deflection angle is then given by. 



Os = Oi - ai{ei) 



(2) 



where 9i is the angular position of the image and 6$ the 
angular position of the source. The amplification matrix at 
any given point is, 
-1 



Defining the generic symmetry matrix. 



(3) 



cos 261 
sin 26 
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- cos 261 



we decompose the amplification matrix as a linear sum: 

= (1 - — SiKp) / — 7c J29c — Si7pjJ2flp. , (4) 

where k is the magnification and 7 the shear. In this frame- 
work, the shear 7 is taken to be a complex number and is 
used to define the quantity 'g as follows: 



Ic + S' 7p. 



1 — Opot 



(5) 



which simplifies in the frame of the perturber j to (neglect 
ing effect of perturber i if i 7^ j): 

7c +T 

9pot\j — 



(6) 



where ffpoJj is the total complex shear induced by the 
smooth cluster potential and the potentials of the per- 
turbers. Restricting our analysis to the weak regime, and 
thereby retaining only the first order ter rns from the lensing 
equation for the shape parameters (see Kneib et al. 1996 ) 
we have: 
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Tl — TS + r pot, 



(7) 



where Tj is the distortion of the image, Ts the intrinsic 
shape of the source, Tpot is the distortion induced by the 
lensing potentials or explicitly in terms of 'g^^t the frame 
of perturber j: 



9i = 9s+gpot\j ^ 9s + 



7c 



+ 



7. 



Pj 



1 K,c l^u 



(8) 



In the local frame of reference of the perturbers, the mean 
value of the quantity 'gj and its dispersion can be com- 
puted in circular annuli (of radius r from the perturber cen- 
ter) strictly in the weak-regime, assuming a constant value 
^^g»9co f-Qj. ^]jg smooth cluster component over the area of 
integration (see Figure 1 for the schematic diagram). 

The result of the integration does depend on the choice 
of coordinate system. In cartesian coordinates (averaging 
out the contribution of the perturbers): 



7c e" 



1 I'^C 



1 — Kc 



7ce 



1 Kc ^Pj 



^2 ""as , 



where 



2iVi 



-I- 



2Ni 



bg 



2Nbg 



(9) 
(10) 

(11) 



o'p(t^) being the width of the intrinsic ellipticity distribution 
of the sources, Nf,g the number of background galaxies av- 
eraged over and cr| the dispersion due to perturber effects 

which should be smaller than the width of intrinsic elliptic- 
ity distribution. In the polar uv coordinates, on averaging 
out the smooth part: 



1 



pj 



7p, 



ss 
2 



+ 



1 Kic f^v 



4 



+ 



■■9p, 



IPj 



1 l^c 



Pj 



where 



p{ts) 



2iV, 



+ 



bg 



2N, 



bg 



(12) 



(13) 



(14) 



From these equations, we clearly see the two effects of the 
contribution of the smooth cluster component: it boosts the 
shear induced by the perturber due to the (kc + i^pj ) term in 
the denominator, which becomes non-negligible in the clus- 
ter center, and it simultaneously dilutes the regular galaxy- 
galaxy lensing signal due to the (t| /2 term (equation 11) 
in the dispersion of the polarization measure. However, one 
can in principle optimize the noise in the polarization by 
'subtracting' the measured cluster signal and averaging it in 
polar coordinates: 



{9 1 -9c)uv 



Ipi 



^Pj 



(15) 



which gives the same mean value as in equation (11) but 
with a reduced dispersion: 



Ss 

2 ■ 



where 



'p(ts) 



2N, 



bg 



(16) 



(17) 



This subtraction of the larger-scale component reduces the 
noise in the polarization measure, by about a factor of two; 
when ~ cr.| , which is the case in cluster cores. Note 
that in subsequent sections of the paper, we plot the aver- 
aged components of r (the quantity measurable from lensing 
observations) computed in the (uv) frame. We reiterate here 
that the calculations above assume that the cluster compo- 
nent is constant over the area of integration (a reasonable 
assumption if we limit our analysis to small radii around the 
centers of perturbers). These results can be easily extended 
to the case when the cluster component is linear (in x and y) 
over the area of integration, the likely case outside the core 
region. This direct averaging prescription for extracting the 
distortions induced by the possible presence of dark halos 
around cluster galaxies, by construction, does not require 
precise knowledge of the center of the cluster potential well. 

2.2 Quantifying the lensing distortion 

To quantify the lensing distortion induced by the individ- 
ual galaxy-scale components using a minimal number of 
parameters to characterize cluster galaxy halos, we model 
the density profile as a linear superposition of two pseudo- 
isothermal elliptical com ponents (PIEMD models derived by 



Kassiola & Kovner 199S) 



E(i?) 



1 



1 



(18) 



with a model core-radius ro and a truncation radius rt 3> 
ro. The useful feature of this model, is the ability to repro- 
duce a large range of mass distributions by varying only the 
ratio rj: defined as ?7 = rt/ro. It also provides the following 
simple relation between the truncation radius and the effec- 
tive radius Rc, rt ^ (4/3)-Ro. Furthermore, this apparently 
circular model can be easily generalized to the elliptical case 
by re-defining the radial coordinate R as follows: 



R^ = 



y 



a + b' 



(19) 



{l + eY (1 

The mass enclosed within radius R for the model is given 
by: 

27rSoro , 



M{R) 



1„ m 



- ro 



and the total mass, which is finite, is: 
Moo — 27rEorort. 
Calculating k, 7 and g, we have, 
ro , 1 



,(20) 



(21) 



ri(R) 



(i-ro/rt) ^(rg + R^) ^(r? + R^ 



--) , (22) 
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Figure 1. Local frame of reference of the perturber: The vector diagram illustrating the choice of coordinate system. The total shear is 
decomposed into a large-scale component due to the smooth cluster and a small-scale one due to the perturbing galaxy. In the frame of 
the perturber, the averaging procedure allows efficient subtraction of the large-scale component as shown in the right panel, enabling the 
extraction of the shear component induced in the background galaxies only by the perturber as shown in the left panel. The background 
galaxies (shown in the left panel of this figure) are assumed to have the same intrinsic ellipticity for simplicity, therefore, we plot only 
the induced components. 



2/to = So 



47rG DuDoi 

C2 Dos ' 



(23) 



where Dis, Dos and Doi arc respectively the lens-source, 
observer-source and observer-lens angular diameter dis- 
tances. To obtain g{R), knowing the magnification k{R), 
we solve Laplace's equation for the projected potential (f>2D, 
evaluate the components of the amplification matrix and 
then proceed to solve directly for f(R), g{R) and t{R). 



ct>2D = 2Ko[^/r^ + R^ - ^ r1 + + (ro - n) In R 

- ro In [r^ + ro y^^^f+RP] + rt In [r? + rt ^ r\ + fl^j] 

(24) 

To first approximation, 



t{R) w 7(i?) = Ko[- 



\ 2_ 

R^ 



+ I^Cv/^^-^o) 



(/R^+7?-n)]. 



Scaling this relation by rt gives for ro < R < rt: 



(25) 



So rt 



a 



where a is the velocity dispersion and for ro < rt < R 



■y{R/rt) oc 



So rt' 



(26) 



(27) 



where Mtot is the total mass. In the limit that R S> rt, we 
have. 



, , 3ko ,2 2i , 2k;o r , 
7(^) = 2R3 1^0 - rt] + ^[n - ro], 



(28) 



and as J? ^ CO, f(R) — » 0, g{R) — » and t{R) — » as 
expected. 



3 RECOVERING GALAXY-SCALE 
PERTURBATIONS 

In this section, we study the influence of the various param- 
eters using the direct averaging procedure on the synthetic 
data obtained fi'om simulations. The numerical simulations 
involve modeling of the global cluster potential, the individ- 
ual perturbing cluster galaxies and calculating their com- 
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by tine smaller 
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potential 
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by the smooth 
cluster 
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> 




X induced by 
the perturber 
potential 

Figure 2. A schematic representation of tlie effect of the cluster on the intrinsic ellipticity distribution of background sources as viewed 
from the two different frames of reference. In the top panel, as viewed in the image frame - the effect of the cluster is to cause a coherent 
displacement x and the presence of perturbers merely adds small-scale noise to the observed ellipticity distribution. In the bottom panel, 
as viewed in the perturbers frame - here the perturber component causes a small displacement t and the cluster component induces the 
additional noise. 



bined lensing effects on a catalog of faint galaxies. We com- 
pute the mapping between the source and image plane and 
hence solve t he lensing equ ation, using the lens tool utility 
developed by Kneib (1993| ), which accounts consistently for 
the displacement and distortion of images both in the strong 
and weak lensing regimes. 



3.1.1 Spatial and Luminosity distribution 

A catalog of cluster galaxies was generated at random with 
the following characteristics. The luminosities were drawn 
from a standard Schechter function with L, — 3.1O^"L0 
and a — —1.25. The positions were assigned consistent with 
the number density !/(r) of a modified Hubble law profile, 



3.1 Modeling the cluster galaxies 



'o 



(29) 
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rt 



(33) 



These imply the foUowing scahng for the Vt/ro ratio -q: 

The total mass M then scales with the luminosity as 
M = 27rEorort 



(34) 



^(cro)^rt = ^cro,Vt,(-^)5+°', (35) 

where a tunes the size of the galaxy halo, and the mass-to- 
light ratio T is given by: 



QO 



-2 




ao 



-2 



T = 12 



(To* 



240 km/ s J \30kpc 



rt* 



-1/2 



(36) 



-2-10 1 
log (R/rJ 



Figure 3. The effect of the assumed scaling relations are exam- 



ined 

g vs. R\rt for various values of (T,/T,*V 0.5, 1.0, 5.0 and 10.0. 



plot of the magnification log 



R/ri and the ohoar 



Therefore, for a = 0.5 the assumed galaxy model has con- 
stant T for each galaxy; if q > 0.5 (q < 0.5) then brighter 
galaxies have a larger (smaller) halos than the fainter ones. 

The physical motivation for exploring these scaling laws 
arises from trying to understand the observed empirical cor- 
relations for early-type (E & SO) galaxies in the fundamental 
plane (FP). The following tight relation between the effec- 
tive radius R^, the central velocity dispersion ctq and the 
mea n surface brightness within Re is found f or cluster galax- 
ies (Jorgensen, Franx, fc Kjacrgaard (1996 ), Djorgovski & 



Davis (1987) and Dressier et al. (1987)) 



The curves on the left panel are for a = 0.5 and on the right 
panel for a = 0.8, (i) solid curves - (L/L*) = 0.5, (ii) dotted 
curves - {L/L*) = 1.0, (iii) short-dashed curves - (L/L*) = 5.0, 
(iv) long-dashed curves - {L/L*) = 10. The magnification is nor- 
malized so that at r = 2 ro, ft = 1; the difference in the slope 
of re abi|ive and buluw log r/n = can be clearly sueu fur bolli 



log Re = 1.24 logcro - 0.82 \og{I)^+cste 



(37) 



that 



can 
epiko appo 



seen 
the plotc of log 



sets of s caling law 
g vs. log {R/rt) at the radius where the mean enclosed surface 
density is approximately equal to the critical surface density. For 
the mass models studied here (cuspy with small core-radii) the 
surface mass density has a large central value and hence a spike 
appears on a scale that is roughly comparable to the core-radius. 



with a core radius ro — 250 fcpc, as well as a more generic 
'core-less' profile of the form: 



"(1-1- 



(30) 



with a scale-radius = 200 fcpc and a = 0.1 which was 
found to be a good-fit to the g alaxy data of the moderat e 
redshift lensing cluster A2218 by Natarajan & Kneib (1996). 
We find however, that the results for the predicted shear 
from the simulations is independent of this choice. 



3.1.2 Scaling laws 

The individual galaxies are then parameterized by the mass 
model of Section 2.2, using in addition, the following scah ngs 
with luminosity (see Brainerd, Blandford, & Small (1996) for 
an analogous treatment) for the central velocity dispersion 
(To, the truncation radius rt and the core radius ro: 



One of the important consequences of this relation is the fact 
that it necessarily implies that the mass-to-light ratio is a 
weak function of the luminosity, ty pically T ~ L" "^ (Jor- 
gensen, Franx, fc Kjaergaard 199(^ ). In terms of our seal 
ing scaling laws, this implies a = 0.8. Henceforth, in this 
analysis we explore both the scaling relations, for a = 0.5; 
the constant mass-to-light ratio case, and a = 0.8; corre- 
sponding to the mass-to-light ratio being proportional to 
L"'^ - consistent with the observed FP. In Figure 3, we plot 
the scaling relations for various values of (L/L*), ranging 
from 0.5 10.0 for a = 0.5 and a = 0.8. Additionally, 
for the constant mass-to-light ratio case, we also plot the 
iso-T curves in terms of the fiducial crj and rj in Figure 4. 
The scaling laws are calibrated by defining an L» (in the 
R-band) elliptical galaxy to have ro» =0.15 kpc, rt, =30.0 
kpc and a fiducial cro,, then chosen to assign the different 
mass-to-light ratios, [ao, = 100, 140, 170, 240, 340, 480 km 
corresponding to T — 2, 4, 6, 12, 24, 48 respectively]. 



3.2 Modeling the background galaxies 

3.2.1 Luminosity distribution 

The magnitude and hence the luminosity for the background 
population was generated consistent with the number count 
distribution measured f rom faint field galaxy s urveys like 
Glazebrook et al. (1994), as well as 



the MDS as reported in 
the more recent results of the number-mag nitude relations 
obtained from the Hubble Deep Field data ( Abraham et al 



(To = (TO: 



ro = ro,( — )■ 



(31) 
(32) 



199C). The slope of the number count distribution used was 
0.33 over the magnitude range m_R = 18 — 26. This power 
law for the number counts implies a surface number den- 
sity that is roughly 90 galaxie s per square arcmi nutes in 
the given magnitude range (see Small et al. (1995b )), which 
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Figure 4. The constant mass-to-light ratio curves are plotted in 
the (aQf,rt-t) plane for an L« galaxy with r) = 200: (i) dot-dashed 
curve - T = 4, (ii) dotted curve - T = 6, (iii) solid curve - T = 
12, (iv) short-dashed curve - T = 24 and (v) long-dashed curve - 
T = 48. 



over the area of the simulation frame [8 arcmin X 8 arcmin] 
corresponds to having ~ 5000 background galaxies. 



3.2.2 Redshift distribution 

The background galaxy population of sources was also gen- 
erated, consistent with the measured redshift, magnitude 
and luminosity distributions (MODEL Z2 b elow) from high- 
redshift s urveys like the A PM and CFRS (Efstathiou et al 



1991 and 



Lilly et al. 1995 respectively). For the normalized 



redshift distribution at a given magnitude m (in the R-band) 
we used the following fiducial forms: 



MODEL Zl. 



N{z,m) = No5{z - 2), 



(38) 



corresponding to the simple case of placing all the sources 
at z —2. 



MODEL Z2. 



simulatior is rriRo =22.0, dZmcdian/rfmfl^O.l a nd Zm pdian^ 
0.58 (see [Brainerd, Blandford, fc Small (199(: ) and Kneit 
et al. (19960; 



3.2.3 Ellipticity distribution 



Analysis o f deep surveys such as the MDS fields (Griffith, 
ct al. 1994 ) shows that the ellipticity distribution of sources 
is a strong function of the sizes of indiv idual galaxies as well 
as their magnitude (Kneib et al. 1996). For the purposes of 
our simulations, since we assume 'perfect seeing', we ignore 
these effects and the ellipticities are assigned in concordance 
with an elli pticity distribution p(rs) deri ved from fits to the 
MDS data ([Ebbels, Kneib, fc Ellis 1997|) of the form. 



p{rs) = rs cxp(-(y) 



= 1.15, 5 = 0.25. 



(41) 



4 ANALYSIS OF THE SIMULATIONS 

We use the above as input distributions to simulate the back- 
ground galaxies and bright cluster galaxies in addition to 
a model for the cluster-scale mass distribution. Analogous 
t o the mass model constructed for the cluster Abell 2218 
( Kneib ct al. (1996 )), we set up an elliptical mass distribu- 



tion for the central clump with a velocity dispersion of 1100 
km s~^ placed at a redshift z = 0.175. The main clump was 
modelled using a PIEMD profile (as in equation (14)) with 
an ellipticity e = 0.3, core radius 70 kpc and a truncation 
radius 700 kpc; therefore the surface mass density of the 
clump falls off as r~^ for r S> »"cut. 

The lens equation was then solved for the specified con- 
figurations of sources and lenses set-up as above and the 
corresponding image frames were generated. The averaged 
components of the shear binned in circular annuli centred 
on the perturbing galaxies was evaluated in their respec- 
tive local (u,v) frames. An important check on the entire 
recovery procedure arises from the fact that by construction 
(choice of the (u,v) coordinate system) the mean value of 
the v-component of the shear < r„ > is required to vanish. 

In the following sub-sections, we explore the dependence 
of the strength of the detected signal on the various input 
parameters. First of all. Figure 5 demonstrates the good 
agreement between the analytic formula for the shear de- 
rived at a given radial distance R produced by a PIEMD 
model as computed in Section 2.2 and the averaged value 
extracted from the simulation on solving the lensing equa- 
tion exactly for the redshift distribution of MODEL Zl. In 
all subsequent plots (Figures 6, 7, 8, 9, 10, 12 and 13) the 
annuli are scaled such that for an L, galaxy, the width of 
each ring corresponds to a physical scale of ~ 20 kpc at 
z = 0.175. 



N{z,m) 



° r(§)zo 



where /? =1.5 and 

Zo = 0.7 Zmcdian H ; (niR - mRo)], 

drUR 



(39) 



(40) 



•Zmodian bciug the median redshift, dz^ncdian/dniR the change 
in median redshift with R magnitude niR. We use for our 



4.1 Error Estimate on the signal 

There are two principal sources of error in the computation 
of the averaged value of the shear aside from the observa- 
tional errors (which are not taken into account in these sim- 
ulations) arising from the effects of seeing etc. (i) shot noise 
(due to a finite number of sources and the intrinsic width 
of their ellipticity distribution) and (ii) in principle the un- 
known source redshifts. Therefore, we require a minimum 
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Figure 5. Demonstrating tiie robustness of the signal extrax;tion 
by comparing the analytic prediction with the measured radi- 
ally averaged shear from the simulation. The signal was extracted 
from a sinmlation run of a PIEMD model with rf» = 30 kpc, rg* 
= 0.15 kpc and velocity dispersion of 480 kms"-*-: solid curve - 
estimate from the analytic formula and overplotted are the mea- 
sured values of the averaged shear. 

threshold number of background objects to obtain a signif- 
icant level of detection. The unknown redshift distribution 
of the sources also introduces noise and affects the retrieval 
of the signal in a systematic way, for instance, the obtained 
absolute value for the total mass estimate for cluster galax- 
ies is an under-estimate for a higher redshift population for 
a given measured value of the shear. The mean (or alterna- 
tively the median) and width of the redshift distribution are 
the important parameters that determine the errors incurred 
in the extraction procedure. 

For the simulation however, we need to obtain an error 
estimate on the signal given that we measure the averaged 
shear for a single realization. In order to do so, the simula- 
tion was set up with a constant mass-to-light ratio (T — 12) 
for the 50 cluster galaxies with 5000 background galaxies, 
and on solving the lens equation the image frame was ob- 
tained. The averaging procedure as outlined in Section 2.1 
was then implemented to extract the output signal with 1000 
independent sets of random scrambled positions for the clus- 
ter galaxies (in addition to the one set of 50 positions that 
was actually used to generate the image); the results are 
plotted as the lower solid curves in Figures 7 & 8. This is a 
secure estimate of the error arising for an individual realiza- 
tion, since this error arises primarily from the dilution of the 
strength of the measured shear due to uncorrelated sources 
and lensed images. We found that the mean error in (tu) 
in the first annulus is 0.040 ± 0.0012 and 0.0048 ± 0.0047 in 




10 20 30 40 50 

r 

Figure 6. Variation of the mean value of the signal in the first 
annulus with mass-to-light ratio T: for T= 2, 4, 6, 12, 24, 48 
of the cluster galaxies plotted for MODEL Zl (solid curve) and 
MODEL Z2 (dotted curve). 

4.2 Variation of the signal with mass-to-light 
ratio of cluster galaxies 

The simulations were performed for mass-to-light ratios (T) 
ranging from 2 ^ 48 (see Figures 6, 7 & 8). The veloc- 
ity dispersion of the fiducial galaxy model was adjusted to 
give the requisite value for T keeping the scaling relations 
intact. The detection is significant for mass-to-light ratios 
T > 4 given the configuration with 50 cluster galaxies and 
5000 background galaxies. The strength of the signal varies 
with the input T of the cluster galaxies, and increases with 
increasing T. As a test run, with T = 0, (i.e. no cluster 
galaxies) and only the large-scale component of the shear, 
we do recover the expected behavior for (tu). The signal 
was extracted for both background source redshift distribu- 
tions MODEL Zl & MODEL Z2. While the amplitude of the 
signal is not very sensitive to the details of the redshift dis- 
tribution of the background population and hence did not 
vary significantly, the error-bars are marginally larger for 
MODEL Z2. This can be understood in terms of the addi- 
tional shot noise induced due to the change in the relative 
number of objects 'available' for lensing; in MODEL Z2 a 
fraction of the galaxies in the low-z tail of the redshift dis- 
tribution end up as foreground objects and are hence not 
lensed, thereby diluting the signal and increasing the size of 
the error-bars marginally. 

4.3 Variation with the number of background 

galaxies 

The efficiency of detection of the signal depends primarily 
on the number of background galaxies averaged over in each 
annulus and therefore on the number that are lensed by the 
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Figure 7. Recovering tiie signal for MODEL Zl: for various val- 
ues of the constant mass-to-light ratio T of the cluster galaxies 
ranging from 2 - 48, (i) lower solid curve - the error estimate (ii) 
upper solid curve - T = 2, (iii) dotted curve - T = 4, (iv) dashed 
curve - T = 6, (v) long-dashed curve - T = 12, (vi) dot-short 
dashed curve - T = 24, (vii) dot-long dashed curve - T = 48. 
Note here that < r„ > is zero as expected by definition of the 
(u,v) coordinate system. 



individual cluster galaxies. For a fixed value of T, the total 
number of background galaxies Nbg was varied, assuming a 
redshift distribution of the form of MODEL Zl. With in- 
creasing Nig, 1000 2500 —^ 5000, the detection is more 
secure and the error does vary roughly as \/ Nbg as shown in 
Figure 9. In principle, the larger the number of background 
sources available for lensing, the more significant the detec- 
tion with tighter error bars: however wo find that a ratio of 
50 cluster galEixies to 2500 background galaxies provides a 
secure detection for T > 4, a larger number of background 
source are required to detect the corresponding signal in- 
duced by lower mass-to-light ratio halos. A secure detection 
in this case refers to the fact that the difference in the mean 
values of the detected signal in the two cases (with Ni,g = 
5000 and Nbg = 2500 background sources) is comparable to 
the mean estimated error per realization computed in Sec- 
tion 4.1. The number count distribution used to generate 
the background sources corresponds to a background surface 
number density of ~ 90 galaxies per square arcmin which 
we find provides a secure detection for T > 4. It is useful 
to point out here that for the standard Bruzual & Chariot 
(95) spectral evolution of stellar population synthesis mod- 
els with solar metallicity, a galaxy that is roughly 10 Gyr old 
(a reasonable age estimate for a galaxy in a a ~ 0.3 clus- 
ter), formed in a single IGyr burst of star formation and 
having evolved passively, one obtains a stellar mass-to-light 
ratio in the R band of ^ 8 with a single power law Salpeter 
IMF with lower mass limit of 0.1 M© and upper mass limit 
125 MQ. With the same ingredients but a Scalo IMF one ob- 



Figure 8. Recovering the signal for MODEL Z2: for various val- 
ues of the constant mass-to-light ratio T of the cluster galaxies 

ranging from 2 - 48, (i) lower solid curve - the error estimate (ii) 
upper solid curve - T = 2, (iii) dotted curve - T = 4, (iv) dashed 
curve - T = 6, (v) long-dashed curve - T = 12, (vi) dot-short 
dashed curve - T = 24, (vii) dot-long dashed curve - T = 48. 



tains a M/L ratio about a factor 2 smaller (~ 4) since there 
arc a smaller proportion of very low-mass stars. Therefore, 
an R-band M/L of 4 for a cluster galaxy is consistent with 
the observed mass just in stars and does not imply the pres- 
ence of any dark mass in the system. Therefore, if dark halos 
were indeed present around the bright cluster members, the 
corresponding inferred mass-to-light ratios would be greater 
than 4, and with 5000 background galax;ies, we would be 
sensitive to the signal as shown in the plots of Figures 6, 7 
& 8. 



4.4 Variation with cluster redshift 

The lensing signal depends on the distance ratio Dis/Dos, 
the angular extent of the lensing objects, the number den- 
sity of faint objects and their redshift distribution. We per- 
formed several runs with the cluster (the lens) placed at 
different redshifts, ranging from z =0.01 to 0.5. We scaled 
all the distances with the appropriate factors corresponding 
to each redshift for both MODELS Zl & Z2. For MODEL Zl 
(Figure 10 and dotted curve in Figure 12), we find that the 
signal (by which we refer to the value of (t„) in the innermost 
annulus) saturates at low redshifts; for 0.01 < iiens < 0.07 
the measurements are consistent with no detection but the 
strength increases as 2iens is placed further away and it re- 
mains significant for upto ziens =0.4, subsequent to which 
it falls sharply once again at 0.5. On the other hand, we 
find that for MODEL Z2 (Figure 11, and solid curve in 
Figure 12), there is a well-defined peak and hence an op- 
timal lens redshift range for extracting the signal. Thus, in 
general, cluster- lenses lying between redshifts 0.1 and 0.3 



10 Priyamvada Natarajan & Jean- Paul Kneib 



0.5 

0.4 
A 0.3 
V 0.2 

0.1 




2 3 
annulus 




2 3 
annulus 




0.3 




Figure 9. Variation of the signal with the number of background 
galaxies for MODEL 7,1: for a given mass-to-light ratio T= 12 of 
the cluster galaxies. We find that the error bars and henec the 
noise decrease as expected with increasing Niig (i) solid curve: 
iVtg = 1000, (ii) dashed curve: N^g = 2500, (iii) dotted curve: 
Nhg = 5000. 



are the most suitable ones for constraining the mean T of 
cluster galaxies via this direct averaging procedure. These 
trends with redsliift can be understood easily, the shear pro- 
duced is proportional to the surface mass density and scales 
as {Dis/Doe) - the saturation at high-redshift is due to the 
combination of two diluting effects (i) the decrease in Di^ 
as the lens is placed at successively higher redshifts (ii) the 
effect of additional noise induced due to a reduction in the 
number of background objects for MODEL Z2. The drop-off 
at low z (for both models) is primarily due to behavior of 
the angular scale factors at low-redshifts. Additionally, the 
shape of these curves is independent of the total mass of 
the cluster (the total mass being dominated by the smooth 
component), therefore even for a subcritical cluster we ob- 
tain the same variation with redshift. 



4.5 Dependence on assumed scaling laws 

In section 3.1, we outlined the simple scaling relations that 
were used to model the cluster galaxies. The choice of the 
exponent a in equation (26) allows the modelling of the 
trends for different galaxy populations: a = 0.5 correspond- 
ing to a constant T and a — 0.8 corresponding to T being 
a weak function of the luminosity. Simulating both these 
cases above, we find that the mean value of the signal does 
depend on the assumed exponent for the scaling law and 
therefore the mass enclosed (Figure 13). We find that while 
the signal is stronger for a — 0.8, since that corresponds 
to the mass-to-light ratio T ~ (;^)°'^ - on average that 
is what we expect compared to the constant mass-to-light 
ratio case; it is not possible however, to distinguish between 



Figure 10. Variation of the signal with cluster redshift for 
MODEL Zl: for a given mass-to-light ratio T= 12 of the clus- 
ter galaxies placing the lens at different redshifts, right panel: (i) 
solid curve: z = 0.1, (ii) dotted curve: z = 0.2, (iii) dashed curve: 
z = 0.3, (iv) long dashed curve: z = 0.4, (v) dot dashed curve: z = 
0.5 and the left panel:(i) solid curve: z = 0.01, (ii) dotted curve: 
z = 0.02, (iii) dashed curve: z = 0.05, (iv) long dashed curve: z = 
0.07, (v) dot dashed curve: z = 0.10 



a correspondingly higher value of the constant M/L and a 
higher value of a. Therefore, the direct averaging procedure 
cannot discriminate between the assumed detail fall-off for 
the mass distribution. 

4.6 Examining the assumption of analysis in the 

weak regime 

While our mathematical formulation outlined in section 2 
is strictly valid only for k <C 1, we examine how crucial 
this assumption is to the implementation of the technique. 
For the output images from the simulations, the magnifi- 
cation K is known at all points. Prior to the averaging, we 
excised the high k regions successively, by removing only the 
lenses in those regions. The results are plotted in Figure 14 
for input T — 12, with the sources distributed as specified 
by MODEL Zl. While the mean peak value of the signal 
does not fluctuate much; on removing the high-K regions, 
we find that the cluster subtraction does get progressively 
more efficient, as evidenced by the sharp fall-off to zero of the 
signal in the second annulus outward. Therefore, while the 
detectability and magnitude of the signal is robust even in 
the 'strong regime', the contribution from the smooth clus- 
ter component, which for our purposes is a contaminant, can 
be 'removed' optimally only in the low k regions. 



5 MAXIMUM-LIKELIHOOD ANALYSIS 
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Figure 11. Variation of the signal with cluster redshift for 
MODEL Z2: for a given constant mass-to-light ratio T= 12 of 
the cluster galaxies placing the lens at different rcdshifts, right 
panel (i) solid curve: z = 0.1, (ii) dotted curve: z = 0.2, (iii) 
dashed curve: z = 0.3, (iv) long dashed curve: z = 0.4, (v) dot 
dashed curve: z = 0.5; left panel (i) solid curve: z = 0.01, (ii) dot- 
ted curve: z = 0.02, (iii) dashed curve: z = 0.05, (iv) long dashed 
curve: z = 0.07, (v) dot dashed curve: z = 0.10 




Figure 12. Variation of the maximum value of the signal with 
redshift: for a given constant mass-to-light ratio T= 12 of the 
cluster galaxies placing the lens at different redshifts for the 2 
background redshift distributions for the sources (i) dotted curve: 
MODEL Zl, (ii) solid curve: MODEL Z2 



Figure 13. Examining the scaling relations - the 2 choices of a 
the exponent of the scaling relation for the truncation radius for 
T= 12. We plot the recovered signal (i) solid curve: a = 0.8, (ii) 
dotted curve: a = 0.5. 




annulus 



Figure 14. The effect of excising the high k regions in the image: 
(for T = 12 of the cluster galaxies) (i) solid curve: k < 0.1, (ii) 
dotted curve:/c < 0.2, (iii) dashed curve:K < 0.3, (iv) long dashed 
curve: K < 0.4, (v) dot dashed curve:K < 0.5. 
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5.1 Limitations of the direct averaging method 

The simulations have enabled us to delineate the role of 
relevant parameters and comprehend the trends with clus- 
ter redshift, the redshift distribution of the sources and the 
mass-to-light ratio of the cluster galaxies. While the direct 
method to estimate the suffers from the following limita- 
tions, specially in the cluster core, (i) being in the 'strong' 
lensing regime, the 'cluster subtraction' is not very effi- 
cient and (ii) the probability of a background galaxy being 
sheared due to the cumulative effect of two or more cluster 
galaxies is enhanced; the core being a region with a high 
number density of cluster galaxies; it does, however, pro- 
vide a robust estimate of the mass-to-light ratio modulo the 
assumed model parameters. 

We now explore applying a maximum-likelihood 
method to obtain significance bounds on fiducial parame- 
ters that characterize a 'typical' galaxy halo in the cluster. 



300 



Schneider & Rix (1996) developed a maximum- likelihood 



prescription for galaxy-galaxy lensing in the field; here we 
develop one to study lensing by galaxy halos embedded in 
the cluster. Schematically, we demonstrate the differences in 
the ellipticity distribution that we are attempting to discern 
in Figure 15. Here we have plotted the intrinsic ellipticity 
distribution of the unlensed sources, sources lensed only by 
a cluster scale component and sources sheared by both a 
cluster scale component and 50 cluster galaxies; from which 
it is obvious that the effect that we intend to measure in 
terms of parameters that characterize the cluster galaxies is 
indeed small, hence recovery of the fiducial parameters in 
this case is considerably harder than in the case of purely 
galaxy-galaxy lensing. 



5.2 Application of the maximum-likelihood 
method 

The basic idea is to maximize a likelihood function of the es- 
timated probability distribution of the source ellipticities for 
a set of model parameters, given the functional form of the 
intrinsic ellipticity distribution measured for faint galaxies. 
We briefly outline the exact procedure below. From the sim- 
ulated image frames we extract the observed ellipticity robs. 
For each 'faint' galaxy j, the source ellipticity can then be 
estimated in the weak regime by just subtracting the lensing 
distortion induced by the smooth cluster and galaxy halos 
given the parameters that characterize both these mass dis- 
tributions, in other words. 



^obs 



Si 7p, 



7c, 



(42) 



where " jp^ is the sum of the shear contribution at a given 
position j from Nc perturbers, and the term 7c is the shear 
induced by the smooth cluster component. In the strong 
regime, similarly, one can compute the source ellipticity us- 
ing the inverse of equation (7). The lensing distortion de- 
pends on the parameters of the smooth cluster potential, 
the perturbers and on the redshift of the observed arclet 
(lensed image), which is in general unknown. Therefore, in 
order to invert equation (7), for each lensed galaxy we need 
to assign a redshift, from a distribution of the form in equa- 
tion (33) given the observed magnitude rrij and take the 
mean of many such realizations. In principle, one needs to 
also correct the observed magnitude for amplification to ob- 




Figure 15. The ellipticity distribution p-rg'- (i) solid curve - in- 
trinsic input ellipticities of the sources, (ii) dotted curve - the el- 
lipticity distribution on being lensed by 50 galaxy-scale mass com- 
ponents and one larger-scale smooth component and (iii) dashed 
curve - the ellipticity distribution induced by lensing only by the 
larger-scale smooth cluster component. 



tain the true magnitude prior to drawing a redshift from 
N{z, m), but this correction in turn depends on the redshift 
as well. An alternative procedure is then to correct for the 
amplification using the median z corresponding to the ob- 
served magnitude from the same distribution. This entire in- 
version procedure is performed within the lens tool utilities, 
which accurately takes into account the non-linearities aris- 
ing in the strong regime. As an input for this calculation, we 
parameterize both the large-scale component and perturb- 
ing galaxies as described in Section 2.2 and Section 3.1 re- 
spectively. Additionally, we assume that a well-determined 
'strong lensing' model for the cluster-scale halo is known. 
For our analysis, we also assume that the functional form 
of p{ts) from the field is known, and is specified by equa- 
tion (34); the likelihood for a guessed model can then be 
expressed as. 



£(ao.,n.,...) =n/°'p(rs^). 



(43) 



However, note that we ought to compute £ for different re- 
alizations of the drawn redshift for individual images (say 
about 10-20) and then compute the mean of the different 
realizations of zj; but it is easily shown to be equivalent to 
constructing the £ for a single realization where the redshift 
Zj of the arclet drawn is the median redshift correspond- 
ing to the observed source magnitude. For the case when we 
perform a Monte-Carlo sum over A^'mc realizations of Zj, the 
likelihood is: 



£(ao.,n„...) =H;'-'H^^'°p(r 



(44) 



where PriTgk) is the probability of the source ellipticity dis- 

j 

tribution at the position j for k drawings for the redshift 
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of the arclet of known magnitude rrij. The mean value for 
A^MC reaUzations gives: 
1 



A'mc 



which written out in integral form is equivalent to 

> j p{Ts^{z))N{z,mj)dz 

(p(rs,)) = 



, m j j dz 



J N{z,'. 



(45) 



(46) 



•^-median 



)) 



«avg being the average redshift corresponding to the magni- 
tude rrij. Therefore the corresponding likelihood £ is then 
simply, 

£ = n,{p{rs^)) (47) 

as before and the log-likelihood / = ln£ = E<^p(Tg^)^. The 
best fitting model parameters are then obtained by max- 
imizing this log-likelihood function I with respect to the 
parameters co* and rt*, the characteristic central velocity 
dispersion and truncation radius respectively. The results of 
the maximization are presented in Figures 16 - 18. For all 
reasonable choices of input parameters we find that the log- 
likelihood function has a well-defined and broad maximum 
(interior to the innermost contour on the plots). The con- 
tour levels are calibrated such that Zmax — I ~ 1,2,3 can 
be directly related to confidence levels of 63%, 86%, 95% re- 
spectively (we plot only the first 10 contours for each of the 
cases in Figures 16 - 18) and the value marked by the dot- 
ted lines denotes the input values. In Figure 16, we plot the 
likelihood contour for the MDS ellipticity distribution (equa- 
tion (34)) - the left panel for an assumed scaling law with 
a = 0.5 and a constant mass-to-light ratio T — 12. On the 
right panel, the corresponding contours for a = 0.8 are plot- 
ted. For the MDS ellipticity distribution, we find that the 
velocity dispersion ctq* can be more stringently constrained 
than the halo size, and the contours are elongated along the 
constant mass-to-light ratio curves and yield an output T 
very nearly equal to the input value. For narrower ellipticity 
distributions both the parameters can be constrained better 
and the inferred T is every nearly equal to the input value. 
We find that there is very little perceptible difference in the 
retrieval of parameters for the two cases with the different 
scaling laws. For a sub-critical cluster (see bottom left panel 
in Figure 18), we find that the parameters are recovered 
just as reliably, which is not surprising and in some sense il- 
lustrates the robustness of the maximum-likelihood method. 
Thus, the physical quantity of interest that can be estimated 
best from the analysis above is the mass M* of a fiducial L* 
galaxy. 

5.3 Estimating the required number of 
background galaxies 

The largest source of noise in our analysis arises due to the 
finite number of objects in the frame. To estimate the re- 
quired signal-to-noise that would permit obtaining reliable 
constraints on both ao* and rt*, we reduced the number of 
background sources to 2500 keeping the number of lenses at 
50 as before. We do not converge to a maximum in the log- 
likelihood, and consequentially, no confidence limits can be 
obtained on the parameters. Therefore, to apply this tech- 
nique to the data we require the ratio r of the number of 



cluster galaxies to the number of background galaxies to be 
roughly r < 0.2, which can be achieved only by stacking 
the data from many clusters. Also, as found from the direct 
averaging procedure, we require ~ 5000 lensed images in or- 
der to securely detect < T > > 4. Although typical HST 
cluster data fields are of the order of [3 arcmin X 3 arcmin] 
have ~ 700 background galaxies (with a 10 orbit exposure) 
of these the shape parameters can be reliably measured only 
for about 200 galaxies, therefore on stacking the data from 
20 (10-orbit) HST cluster fields, we shall be able to constrain 
statistically the mean mass-to-light ratios as well as the 2 
fiducial parameters. 



5.4 Uncertainties in the smooth cluster 
component 

In all of the above, we have assumed that the parameters 
that characterize the smooth cluster-scale component are 
very accurately known which is unlikely to be the case for 
the real data. We investigate the error incurred in retrieving 
the correct input parameters from not knowing this central 
strong lensing model well enough. So we can now place limits 
on the order of magnitude of errors that can be tolerated 
due to the lack of knowledge of the exact position of the 
cluster center and the velocity dispersion of the main clump. 
In Figure 18, we see that an uncertainty of the order of 
20 arcsec in the position of the center yields unacceptably 
incorrect values for tJo, and ru- Conversely, if the center 
is off by only 5 arcsec or so, for both the critical cluster 
and the sub-critical one, the results remain unaffected and 
we obtain as good a retrieval of the input as when the 
position of the center is known exactly. Similarly, in Figure 
17, we demonstrate that an error of ~ 5% in the velocity 
dispersion is enough to make the max-likelihood analysis 
inconclusive, but an error of ~ 2-3% at most would still 
enable getting sensible bounds on both parameters. 



6 CONCLUSIONS AND PROSPECTS 

We conclude this section and assert that both the maximum- 
likelihood method and the direct averaging method devel- 
oped in this paper can be feasibly applied to the real data 
on stacking a minimum of 20 WFPC2 deep cluster fields. 
These methods are well-suited to being used simultaneously 
as they are somewhat complementary, both yield the sta- 
tistical mass-to-light ratio reliably and while the averaging 
does not require either the knowledge of the center or any 
details of the strong lensing model, it also cannot provide 
the decoupling of the 2 fiducial parameters and hence no 
independent constraints on the velocity dispersion and the 
halo size can be obtained, meanwhile the maximum likeli- 
hood approach permits estimation of the fiducial ao* and rt* 
{ao* more reliably than rt*), it necessarily requires knowl- 
edge of the cluster center and the central velocity dispersion 
rather accurately. In offset fields, however, where the gra- 
dient of the smooth cluster potential is constant over the 
smaller scales that we are probing, we expect both methods 
to perform rather well. 

In this paper, we have not investigated the likely sources 
of error in the real data, which we do in detail in a subse- 



quent paper (Natarajan, Kneib, fc Small 1996), our Simula- 
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Figure 16. Log-likelihood contours for the retrieval of the fiducial parameters cro* and rt* - the input values are indicated by the 
intersection of the dashed lines. In the left panel: For the MDS ellipticity distribution, with assumed scaling a = 0.5, right panel: the 
same with a = 0.8 



tions have enabled the study of the feasibility of application 
to HST cluster data in as far as a statistical estimate of 
the required number of background galaxies required for a 

significant detection given the limitations in the accuracy 
to which the input parameters (like the strong lensing mass 
model and hence the magnification) are presently known. 
Our analysis points to the fact that the extraction of the 
signal would therefore be feasible if approximately 20 - 25 
clusters are stacked, and the enterprise is specially suited 
to using the new ACS (advanced camera for survey) due to 
be installed on HST in 1999. Additionally, since there ex- 



ists a well-defined optimum lens redshift for signal detection 
(0.1 < 2iens <0.3), it might be useful to target clusters in 
this redshift range in future surveys in order to apply the 
techniques developed here. In our proposed analysis with 
the currently available HST data, we intend to incorporate 
parameters characterizing the smooth cluster (main clump) 
alongwith those of the perturbing galaxies into the maxi- 
mum likelihood machinery. 

In summary, we have presented a new approach to in- 
fer the possible existence of dark halos around individual 
bright galaxies in clusters by extracting their local lensing 
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Figure 17. Sensitivity of log-likelihood contours to the strong lensing input parameters: examining the tolerance of the significance 
bounds obtained on (Tq* ^■iid ^t* with regard to the accuracy with which the cluster velocity dispersion needs to be known. All plots 

arc for input T = 12, a = 0.5 and the MDS source ellipticity distribution. Top left panel: given that the exact value of the velocity 
dispersion is known (value in this case is 1090 kms~l), top right panel: the velocity dispersion known to within 2%, bottom left panel: 
attempt to retrieve the incorrect scaling law - input a = 0.5, log-likelihood maximized for a = 0.8, bottom right panel: retrieval with 
fewer background galaxies 



signal. The composite lensing effect of a cluster is modeled in 

numerical simulations via a large-scale smooth mass compo- 
nent with additional galaxy-scale masses as perturbers. The 
correct choice of coordinate frame i.e. the local frame of each 
perturber, enables efficient subtraction of the shear induced 
by the larger scale component, yielding the averaged shear 
field induced by the smaller-scale mass component. Cluster 
galaxy halos are modeled using simple scaling relations and 



the background high redshift population is modeled in con- 
sonance with observations from redshift surveys. For several 
configurations of the sources and lens, the lensing equation 
was solved to obtain the resultant images. Not surprisingly, 
we find that the strength of the signal varies most strongly 
with the mass-to-liglit ratio of the cluster galaxies and is 
only marginally sensitive to the assumed details of the pre- 
cise fall-off of the mass profile. We also find that there is 
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Figure 18. Sensitivity of the log-likelihood contours to input parameters: examining the tolerance of the significance bounds obtained 
on (TO* and rt» given the accuracy to which the cluster center needs to be known. All plots are for input T = 12, a = 0.5 and the MDS 
source ellipticity distribution. Top left panel: knowing the cluster center exactly for the critical cluster, top right panel: knowing the 
center to within 5 arcsec, bottom left panel: knowing the center exactly for the sub-critical cluster and the bottom right panel: for the 
subcritical cluster, center known to within 5 arcsec. 



an opti mum lens redshift range for detection of the signal. 



Althouj;h the entire procedure works in the 'strong lensing' 
regime as well, it is less noisy in the 'weak regime'. The 
proposed maximum-likelihood method independently con- 
strains the halo size and mass of a fiducial cluster galaxy 
and we find that the velocity dispersion and hence the mass 
of a fiducial galaxy can be more reliably constrained than 
the characteristic halo size. Examining the feasibility of ap- 
plication to real data, we find that stacking ^ 20 clusters 



allows a fir st attempt at extraction (Natarajan, Kneib, & 



Smail 1996). The prospects for the application of this tech- 



nique are potentially promising, specially with sufficient and 
high-quality data (either HST images or ground-based ob- 
servations under excellent seeing conditions of wider fields); 
the mass-to- light ratios of the different morphological/color 
types in clusters for instance can be probed. More impor- 
tantly, comparing with similar estimates in fields offset from 
the cluster center would allow us to make the essential con- 



Lensing by galaxy halos in clusters of galaxies 17 



nections in order to understand the dynamical evolution of Tyson J. A., Valdcs F., Jarvis J. F., Mills A. P., 1984, ApJL, 281, 
galaxies in clusters and the re-distribution of dark matter 59 

within smaller scales within clusters. Application of this ap- West M., Bothun G., 1990, ApJ, 350, 36 
proach affords the probing of the structure of cluster galaxies 
as well as the efficiency of violent dynamical processes like 
tidal stripping, mergers and interactions which modify them 
and constitute the processes by which clusters assemble. 
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